Nonequilibrium spintronic transport through an artificial Kondo impurity: 
Conductance, magnetoresistance and shot noise 
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We investigate the nonequilibrum transport properties of a quantum dot when spin flip processes 
compete with the formation of a Kondo resonance in the presence of ferromagnetic leads. Based upon 
the Anderson Hamiltonian in the strongly interacting limit, we predict a splitting of the differential 
conductance when the spin flip scattering amplitude is of the order of the Kondo temperature. We 
discuss how the relative orientation of the lead magnetizations strongly influences the electronic 
current and the shot noise in a nontrivial way. Furthermore, we find that the zero-bias tunneling 
magnetoresistance becomes negative with increasing spin flip scattering amplitude. 

PACS numbers: 72.15.Qm, 72.25.Mk, 73.63.Kv 
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Introduction. — The Kondo effect in quantum dots 
(QD's) occurs because of a strong antiferromagnetic cou- 
pling between the conduction band electrons in the elec- 
trodes and the localized spin in the QD through higher- 
order tunneling processes. The resulting correlated mo- 
tion gives rise to a Kondo singularity in the quasiparticle 
density of states at the Fermi level Ep Q ■ As a conse- 
quence, the conductance is enhanced below the Kondo 
temperature Tk 2]. The fact that the parameters that 
define Tk are fully controllable and the ability to ap- 
ply external fields in QD's have spurred a good deal of 
works addressing nonequilibrium situations 0. At the 
same time, recent advances in the control and manip- 
ulation of the spin degree of freedom in magnetic en- 
vironments 0| has paved the way for startling poten- 
tialities in spintronic devices (fj and quantum computa- 
tion jfj. A single spin- 1/2 confined in a discrete level 
reaches the ultimate degree of miniaturization in semi- 
conductor physics. Thus, it is crucial to understand the 
properties of spintronic transport in QD's. The main 
focus has so far been to investigate magnetic effects in 
QD's in the Coulomb blockade regime UQ- In compar- 
ison, the issue of spin-dependent physical effects in the 
Kondo resonance, which seems to be much richer, has 
been relatively unexplored. In particular, one might ask 
to what extent the Kondo cloud is sensitive to the lead 
magnetization. Or, is the many-body Kondo state robust 
against spin flip processes? Answering these questions is 
not only interesting for QD's but also for similar coherent 
phenomena in molecular transistors jjj. 

Consider first a noninteracting QD with a single en- 
ergy level £o attached to two ferromagnetic contacts. In- 
tradot spin flip processes lift the level degeneracy, yield- 
ing e ± R, where R is a phenomenological spin flip scat- 
tering amplitude. When R is greater than the tunnel- 
ing induced broadening T, the linear conductance will 
show a splitting sensitive to the lead polarization align- 
ment. Similarly, in the strongly interacting case we ex- 
pect the Kondo induced spin fluctuations to become af- 
fected by the combined influence of ferromagnetic leads 



and spin flip processes. Our main findings reflect this 
physics, summarized in Figs. ^ and [2I a ) : (i) Let Tk 
denote the Kondo temperature of a QD for a given 
magnetic configuration (i? ^ and lead magnetization 
7] 7^ 0). For R/Tk < 1 the differential conductance 
Q = dl/dVdc (Kic is the bias voltage) shows the dis- 
tinctive zero bias anomaly (ZBA) of Kondo physics 
whereas for R/Tk > 1 the ZBA smoothly splits at finite 
Vdc into two peaks separated by a distance 8. With de- 
creasing Tk (e.g., by lowering e with a gate voltage), the 
spin flip scattering is shown to suppress the Kondo state, 
(ii) We find that the exact form of Q depends on whether 
the relative orientation of the lead magnetizations is par- 
allel (P) or antiparallel (AP). In addition, in the P case 
the ZBA splitting develops at different R with regard to 
the AP case. We find that in the AP configuration the 
variation of the Kondo temperature with R is indepen- 
dent of r\. Strong effects of R and 77 in the measurement 
of the magnetoresistance and the current fluctuations at 
T = (shot noise) are predicted as well. 

Model. — There exist in the literature very scarce theo- 
retical investigations that deal with Kondo transport in 
QD's in the presence of ferromagnetic leads jlfi IllL Il2j. 
To the best of our knowledge, no complete picture of the 
Fermi- liquid behaviour (i.e., at temperatures T -C Tk) 
of the Kondo effect in ultrasmall magnetic tunnel junc- 
tions has been yet put forward. We consider a QD (re- 
gion 0) with large on-site Coulomb interaction U cou- 
pled via Vk a to two Fermi-liquid reservoirs labeled with 
a = {L, R} with chemical potentials [II and We 
describe the system with a N = 2 fold degenerated An- 
derson Hamiltonian allowing for intradot spin flips. In 
the limit U — > 00, QD double occupancy is forbidden 
and the slave-boson representation may be used to 
write the Hamiltonian as follows: 

W = £w4aa C W + £ of°f° + ( R f\fi + H - c 
kacr <y 

+-7^ EWk"<w bt /a + h.c.) + A(&t 6 + j2 flu - m 

kaa o 
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where c\ a(J (cfc QCT ) is the creation (annihilation) operator 
for an electron in the state k with spin a = {t:4} in 
the lead a. Ferromagnetism in the leads arises through 
spin-dependent densities of states v a . We have assumed 
implicitly that both ferromagnetic leads possess the same 
(z) easy axis. The generalization to noncollincar mag- 
netic moments is straightforward [ljj. I n Eq. the 
spin flip term is assumed to be coherent in the sense 
that R does not involve spin relaxation since each flip- 
flop process can be reversible ^ij. Only when r\ = (no 
magnetization) do we include in the QD a vanishingly 
small Zeeman splitting to intentionally break the SU(2) 
spin symmetry (e 0T - £ i = A z -> 0+). Finally, fl (/ CT ) 
is a pseudofermion operator that creates (annihilates) a 
singly occupied state and the auxiliary boson operator 
(b) creates (annihilates) an empty state in the QD. The 
last term is a constraint enforced by the replacement of 
the QD second-quantization operators by the /'s and 6's 
when U — > oo with an associated Lagrange multiplier 
A 0. 

The solution of the Hamiltonian 0) can be found in 
the mean field approach, which is the leading order in 
a 1/N expansion, extended to deal with nonequilibrium 
situations 0|. This way one sets b(t)/y/N to a c-number 
corresponding to its expectation value b = (b(t))/\/~N, 
thereby neglecting charge fluctuations. This is correct as 
long as we are interested only in spin fluctuations. Now, 
by calculating the equation of motion for b and taking 
into account the QD constraint we arrive at: 

^2v ka G< afia (t,t) = -iNX\b\ 2 , (2a) 
J2G< afia (t,t) = t(l-N\b\ 2 ), (2b) 

a 

which self-consistently determine the unknowns \b\ 2 and 
A. We have defined Vka = |&| 2 Vfc Q . The nondiago- 
nal (in the layer indices) Green function G< a0a (t,t) = 

*( C L<7 (*)/*(*)) Can be CaSt in terms 01 ^Oo-.OctC*)*) = 

i(fl(t)f a (t)), with the help of the equation of motion 
of the operators and then applying the analytical con- 
tinuation rules in a complex time contour |17| . After 
lengthy algebra, the Fourier transform of 0(T becomes 
G< ai0a (s) = 2iZ a Ue)(m\ 2 fl + fff%)/\MSM§- 
R 2 \ 2 , where f a (e) is the Fermi function of lead a and 
Mj?° = e — ioa ± s^ia ^a- Notice that by virtue of 
Kondo correlations, the original energy level Eq g is trans- 
formed into £o<r = £o<r + A- Likewise, the linewidths 
Y a a = |6| 2 T^ renormalize the bare couplings Y a a (e) = 
(tu/2) J2 ka \Vka\ 2 $(£-£ka<j)- In what follows, T a a is taken 
as T£ = T a a (E F ) for -D < e < D (D is the energy cut- 
off). Within the slave-boson mean-field theory, Kondo 
physics arises when the auxiliary boson field b is con- 
densed. At that point, gives roughly Tk and A shifts 
the resonant level up to Ep. 

The magnetization at reservoir a is r\ a — (v\ — 




FIG. 1: (a) Differential conductance Q versus bias voltage Vdc 
in the parallel case with lead magnetizations rjL = T]R = f] = 
0.5 for four different values of the spin flip amplitude 7?. (b) 
Dependence of the Q—Vdc curves on r\ for R — 0.75T^. Inset: 
Kondo temperature of the QD (Tk ) as a function of R for the 
lead polarizations depicted in (b). The crossing points with 
the straight line Tk = R marks the transition to a nonzero 
splitting in Q(Vd c ). 

+ As -1 < 7? < 1, the linewidths f£ 

will be generally spin dependent: = (1 ± r} a )T a /2. 
Hereafter we shall deal with symmetric leads so that 
Tl = = T/2. Notice that spin flip effects are fully 
included due to the presence of R in the denominator 
of GQ a0a (e). Hence, the closed form of Gq g 0a {e) given 
above represents the distribution function of a QD in the 
presence of spin flips, ferromagnetic leads and Zeeman 
splitting . For R = and rj a = we get at equilibrium 
the position of the Kondo resonance of the impurity at 
£o = £o| = £oi with width r a = L^ + 1^ and Kondo tem- 
perature T£ = (el + T 2 ) 1 / 2 = Dexp(-7r|e |/r). From 
the solution of Eq. JSJl, we can obtain the current travers- 
ing the dot: / = (-e/wh)hn^ a fj J[2G r 0afia (e)f L (e) + 
Gq ct Qa (e)]de, f(E) being the Fermi function. In the fol- 
lowing, we shall present mainly nonequilibrium results 
for fj, L = — fi R = eVd c /2, £ = — 3.5r and D — —60T 
(Tpr w 10 _3 r). The reference energy is set at Ep = 0. 

Results. — Let us first focus on the case of polarized 
reservoirs with parallel (P) alignment (rji, — tjr = v)- 
Figure ^a) shows Q for partial lead magnetization (r\ = 
0.5) and different values of R. As seen, Q is a direct mea- 
sure of how the spin flip mechanisms weaken the Kondo 
effect. For relatively small values of R (R < Tk) the 
height of the ZBA decreases but Q does not split. The 
appearance of a splitting (5^0 occurs at R ~ Tk- Im- 
portantly, due to that coexistence of the Kondo state 
and the intradot spin flips, the positions of the peaks 
in Q are renormalized by Kondo correlations, i.e., these 
peak positions are not trivially related with their single- 
particle counterparts (which should be centered around 
Vdc = ±2i?). For R > T K , the mean field approach pre- 
dicts a complete quenching of the ZBA (Tk = 0) because 
the solution of Eq. @ is the trivial one (6 = 0), which 
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FIG. 2: (a) Q-V^c curves in the antiparallel case for r)L = 
— TjR — 0.5 and different values of R. Inset: Tk as a function 
of R. Notice that this curve is independent of rj so that the 
transition to the splitting in Q(Vdc) takes place at the same 
point independently on the value of rj. (b) Tunneling magne- 
toresistance M at rj — 0.5 as a function of the bias voltage 
Vdc for the first four values of R plotted in Fig.EJa). 



is unphysical. A more detailed theory including charge 
fluctuations (noncrossing approximation) would correct 
this. They price that one has to pay, however, is the 
appearance of spurious peaks in the conductance. Only 
for illustrative purposes have we also plotted in Fig. Of a) 
the case with R = 1.25T& 0. Here, Q is maximal at 
Vd c ~ ±2i2. Notice that for R > T K our model overes- 
timates the height of the peaks as they should quickly 
tend to zero in the Coulomb blockade valley. 

In Fig. IHb) we plot the dependence of Q on the lead 
polarization (for P alignment) at a fixed value of R. The 
overall effect of 77 is to diminish Q. In particular, the 
ZBA quenching can be ascribed to a weaker Kondo state 
because of a smaller T%. The Kondo temperature can 
be easily obtained from Eq. |l3| but for R ^ and 
rj 7^ we find an implicit equation that can be worked 
out numerically (in contrast to the AP case, see below). 
We depict the results in the inset of Fig.^b). The split- 
ting 5^0 takes place at a lower R upon increasing 77, in 
agreement with Fig.^b). We have checked that our nu- 
merical results fulfill two requirements, namely: (i) they 
are independent of the sign of 77, for Tk must be invariant 
after flipping simultaneously all the electron spins in both 
reservoirs, and (ii) with increasing 77, Tk decreases [see 
inset in Fig. Q^b) along the vertical axis R = 0] and even- 
tually vanishes for full lead polarization (77 — > 1), i.e., no 
Kondo effect may arise in the case of half-metallic leads. 

Unlike the P case, for the AP alignment (t]l = —rjn, = 
rj) the ZBA splitting is generated at larger values of R 
[see Fig. |3fa) for 77 = 0.5]. Nevertheless, we observe that 
the ZBA height is reduced even for R = 0. The reason 
is that Q decreases roughly by a factor 1 — rj 2 when the 
magnetic moments point along the opposite directions. 
More interesting is the fact that the ZBA width is given 
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FIG. 3: Fano factor versus bias voltage as a function of R 
for the parallel alignment with rjL ~ rjR = 0.5. (b) Voltage 
dependence of the shot noise for the unpolarized, parallel and 
antiparallel case at R = 0. Inset: transmission probability as 
a function of energy (the origin is taken at Ef = 0) for the 
same values of 7? as in (a). 



by a 1 77 1 -independent Kondo temperature. We find Tk — 
V (Tk) 2 ~ i- e -j Tk does not depend on [77! [see inset 
of Fig. Efa)]- The underlying physics for this is that 
both spin channels are coupled to the QD in the same 
way (i.e., t\ + = T L L + T l L = T) so that T K persists 
unaltered regardless of the value of 77. As a result, the 
splitting in Q(Vd c ) takes place always at the same point 
independently on the degree of lead magnetization. 

We have found as well interesting features in the tun- 
neling magnetoresistance M = (Q p — Q AP )/Q AP [see 
Fig. Efb)]. For definiteness, let us focus on the Vdc = 
case. In conventional magnetic tunnel junctions, it is 
usual that the P current is larger than the AP current, 
thereby resulting in an increase of the impedance of the 
system when the co nfig uration is switched from the P to 
the AP alignment 20] . This is really the case for the 
Kondo ZBA in the absence of spin flips [see Fig.[2Ib) for 
R = 0]. However, increasing R does not lead only to a 
reduction of M but also to a reversal of the sign of M. 
This arises from the fact that the Kondo effect is more 
robust against spin flips in the AP configuration than in 
the P alignment, as stated above. 

Shot noise. — As known, current fluctuations due to 
charge granularity (shot noise) can provide informa- 
tion additional to the averaged current [21]. The shot 
noise is given by the current-current correlation func- 
tion S{t - t') = (Si(t)SI{t')), where §1 = 1- (I) mea- 
sures the current fluctuations. (Because we consider 
a two-terminal system, we have dropped the lead in- 
dices.) An important figure of merit is the Fano factor: 
7 = P(0)/2e(7), where P(0) is the noise power spectrum 
at zero- frequency [i.e., the Fourier transform of S(t — t')]. 
The Fano factor determines the deviations of the shot 
noise away from the Poissonian value of a classical con- 
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ductor due to, e.g., strong electron-electron interactions 
as those giving rise to the Kondo effect. Figure shows 
the influence of spin flips (a) and lead magnetization (b) 
in 7. At low bias, 7 behaves as f — T(Ep) [2lj . where T 
is the transmission probability. As shown in Fig. [3f a) for 
R = 0, the Kondo unitary limit is reached for a partial 
lead magnetization in the parallel case. Because of the 
correlated motion of the electrons which lead to the sin- 
glet formation there is a suppression of 7 at zero bias. As 
R is turned on, the transmission departs from its unitary 
limit [see inset of Fig. EJU)] because the Kondo effect is 
quenched and 7 takes on a nonzero value. Finally, for the 
largest R the Fano factor tends to its Poissonian limit 
as spin flip processes reduce the Kondo resonance and 
a two-peak structure arise in the transmission probabil- 
ity. Thus, while the splitting of Q requires large applied 
voltages to be observable, the shot noise shows the same 
effect already in the zero voltage limit. On the other 
hand, the behavior of 7 at larger Vd c depends on the 
specific transport mechanism. At R = 0, correlations be- 
tween the band electrons and the localized electron dom- 
inate so that 7 increases with Vd c [12] • For R/T 1 ^ — > 1, 
the Kondo resonance is split [see inset in Fig. C^b)] and 
transport is maximized as Vdc sweeps across the new res- 
onances. For comparison between normal metals and fer- 
romagnetic electrodes, we observe from Fig.^Jb) (R = 0) 
that (i) 7 increases more rapidly for the P than for the 
unpolarized case due to the 77 dependence of Tr-, and (ii) 
for AP alignment the ZBA diminishes by a factor f — ?y 2 
and therefore we find 7 ~ rj 2 at low bias. 

Conclusion. — Using a slave-boson mean-field theory, 
we have shown that the combined influence of ferromag- 
netic electrodes and spin flip transitions in the Kondo 
physics of a QD manifests itself in the nonequilibrium 
transport properties of the system. The Kondo temper- 
ature (which can be regarded as the coupling strength of 
the singlet formation) is shown to be suppressed with in- 
creasing lead magnetization in the parallel alignment but 
the Kondo state is remarkably stable in the antiparallel 
case. This may lead to negative tunneling magnetoresis- 
tance at zero-bias when the spin flip scattering amplitude 
is enhanced. The overall behavior is confirmed with shot 
noise calculations. Experimentally, we believe that the 
effects addressed in this paper should be visible within 
the scope of present techniques as our energies are within 
the Kondo scale. 
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